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Abstract 

We study the impacts of nonuniversal Z' model, providing flavor changing neutral current at 
tree level, on the branching ratios (BRs), CP asymmetries (CPAs) and polarization fractions of 
B decays. We find that for satisfying the current data, the new left- and right-handed couplings 
have to be included at the same time. The new introduced effective interactions not only could 
effectively explain the puzzle of small longitudinal polarization in i? — > K*(j) decays, but also provide 
a solution to the small CPA of n^K^. We also find that the favorable CPA of B^ tt^K^ 

is opposite in sign to the standard model; meanwhile, the CPA of B^ tt^K has to be smaller than 
—10%. In addition, by using the values of parameters which are constrained by i? ^ irK, we find 
that the favorable ranges of BRs, CPAs, longitudinal polarizations, and perpendicular transverse 
polarizations for {B^ p^K*, Ba p^K*^) are (17.1 ±3.9, 10.0 ±2.0) x 10-^ (3 ±5, 21 ±7)%, 
(0.66 ± 0.10, 0.44 ± 0.08) and (0.14 ± 0.10, 0.25 ± 0.09), respectively. 



* Email: physchen@mail.ncku.edu.tw 
Email:hatanaka@phys.nthu.edu.tw 



1 



I. INTRODUCTION 



Some puzzles have been found in B meson decaying processes, such as (a) the large 



(BRs) of 5 ^ Ki (b) the small longitudinal polarizations oi B ^ 

and (c) the unmatched CP asymmetries (CPAs) and BRs in B ttK 



branching ratios 
K*(f) decays |3l 14^ 

decays 0, Isl. S llOL lllj . The interesting thing is that the mazy problems are all related to 
the flavor changing neutral current (FCNC) b ^ s decays. The resultants push us not only 
to consider more precise QCD effects but also to speculate the existence of new physics. 

It is known that the FCNC processes usually arise from the loop corrections in the 
standard model (SM) like models. The preference of loop corrections originates from the 
strict constraints of K meson oscillation. However, the constraints of K system are only on 
the first two generations, the FCNC at tree level in third generation has no significant limit 
yet. In especial, the direct constraint from Bg — Bg mixing, dictated by 6 — s interactions, 
is a lower bound in experiment 0]. Hence, it will be interesting to investigate the models 
which FCNC occurs at tree level and they could provide the solutions to the puzzles in 6 ^ s 
decays. 

One of simple extensions of SM for the effects of FCNC is the unconventional Z' model, 
in which FCNC is arisen from the family nonuniversal couplings i.e. the couplings of Z' 
to different families are not the same. One of possible ways to get the family nonuniversal 
couplings is to include an additional U{iy gauge symmetry; and then by the requirement 
of anomaly free, the gauge charges for different families are different |l2| . Other models 
giving the family nonuniversal Z' interactions could be referred to Ref. jl3[. The detailed 
phenomenological analyses for various low energy physics could be found in Ref. Es- 
pecially, the implications on time-dependent CPA of i? ^ Ks4> and on the BRs of B —>■ rj'K 
have been studied by Ref. ^|. Moreover, the solution to the B — > ttK puzzle by the nonuni- 
versal Z' couplings is also discussed by the authors of Ref. |16i] . To further pursue the effects 
on B decaying processes, in this paper, we will take all the measurements of -B ^ K^*^(j) and 
B — > TiK into account to constrain the free parameters. By the constrained parameters, we 
investigate the implications of the Z' model on BRs, CPAs and polarization fractions (PFs) 
for the decays B K*(f), B — * ttK and B pK*. 

For two-body color-allowed processes of B decays, it is known that the dominant hadronic 
effects are the factorized parts which could be simply expressed as the multiplication of 
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effective Wilson coefficients (WCs), the decay constant, and the transition form factors, 
e.g. A oc C{fi)fMiF^^'^^{q'^ = m|^J. Furthermore, by the observed CPA of 0(10%) in 
Bd 71^ K^, we know that the large strong phases have to be introduced. In order to 
self-consistently calculate the hadronic effects, we employ the perturbative QCD (PQCD) 
approach jl?! . isl to evaluate the hadronic matrix elements, where the large strong phases 
could be generated by the annihilation effects of the effective operator {V — A) ^ {V + A). 

It is known that in B system which is composed of a heavy quark and a light quark, the 
residual momentum of the light quark is typically k ~ 0{mB ~ frih) with 'mB{b) being the 
mass of the B-meson (b-quark). In B decays, if we regard that the processes are dominated 
by short-distant interactions, for catching up the energetic quark from the b-quark decay to 
form a energetic meson with the typical energy being ©(m^), the light quark inside the B 
meson has to obtain a large momentum from the b-quark via the gluon exchange. Hence, 
the momentum transfer carried by the hard gluon could be estimated to be ki — k2, where 
ki and k2 denote the momenta of spectator quarks inside the B meson and produced meson, 
respectively. In terms of light-cone coordinates, the large components of ki (z = 1,2) could 
be defined by k^ = Xim^/v^ with Xi being the momentum fractions. Hence, the squared 
momentum of the exchanged hard gluon is = XiX2Tn\. As known that the residual 
momentum of light quark in the B meson is 0{mB — fn^), Xi is roughly 0{mB — m}^)/mB. 
Since the produced light meson is energetic, the momenta of valence quarks should be 
0{mB/\/2), i.e. X2 ~ 0(1). By taking xi = 0.16, X2 = 0.5 and niB{b) = 5.28(4.4) GeV, we 
get ^/cp' ~ 1.5 GeV. Since the value reflects the typical reacting scale of B decays in the 
framework of the PQCD, for the SM contributions, in our calculations the values of weak 
WCs are estimated at the scale /i ~ 1.5 GeV. 

The paper is organized as follows: In Sec. |Hj we introduce the nonuniversal Z' effects 
for h ^ s transition. In Sec. IIIH based on the flavor diagrams, we explicitly write out 
the factorizable amplitudes associated with the new physics for the decays Bd K^*^(j), 
B ttK and B pK*. In addition, we also define direct CPA and PFs. Then by setting 
the values of parameters, in Sec. lIVI we give the calculated values for hadronic effects, present 
various current experimental data for constraining the unknown parameters, display the SM 
predictions and discuss the results of the Z' model. Finally, we give a summary. 



II. FCNC FOR s TRANSITION IN THE Z' MODEL 

In this section, we will introduce the neutral current interactions in the SM and its 
extension with an extra Z' boson. Since we will study the nonleptonic decays, in following 
discussions we only concentrate on the quark sector. Although we concentrate on the study 
of new physics, the used notation for new interacting operators will be similar to those 
presented in the SM. Therefore, it is useful to introduce the effective operators of the SM. 
Thus, we describe the effective Hamiltonian for b sqq decays as 
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q=u,c 



1=3 



(1) 



where Vg = V*gVgb are the Cabibbo-Kobayashi-Maskawa (CKM) [20j matrix elements and 
the operators Oi-Oio are defined as 

Oi^ = {Saqi3)v-A{qpK)v~A , O^''^ = {Saqa)v-A{qi3bi3)v-A , 

O3 = iSaba)v-A ^{qf3qf3)v-A , O4 = {Sab/3)v-A ^(^/jgajv-A , 

O5 = {Saba)v-A ^^(^/^g/?) V+A , = {Sab,3)v-A ^{qi3qa)v+A , 

q g 

Or = ^(sa&a)y-A eq{qf3qf3)v+A , Og = -{sabp)v~A eg{qf}qa)v+A , 



O9 = :^{Saba)v-A^eg{qpq(^)v-A , Oio = -{Sabf3)v-A^eq{qfsqa)v-A 



(2) 



g g 
with a and (3 being the color indices. In Eq. (^, O1-O2 are from the tree level of weak 
interactions, O^-Oq are the so-called gluon penguin operators and Oj-Oiq are the electroweak 
penguin operators, while Ci-Cio are the corresponding WCs. Using the unitarity condition, 
the CKM matrix elements for the penguin operators O3-O10 can also be expressed as \4 + 
Vc = -Vt. 

For studying the Z' model, as usual we describe the Lagrangian for the neutral current 



interactions in terms of weak eigenstates as 



J'2 



7L 



qLJij-^L 



Pl + e 



~R\.. 



(3) 
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where the subscript of qi denotes flavor index of quark, Pl^r = (1 =F 75)/2, and are 
the neutral gauge bosons, corresponding to SU{2) x f/(l) and the extended Abelian gauge 
symmetries, and gi and g2 are the associated gauge couplings, respectively. We note that 
^L,R universal couplings of the SM while e^'^ are 3x3 matrices and denote the effects 
of nonuniversal couplings. In general, the and Z2 bosons will mix each other so that the 
physical states of Z bosons could be parametrized by 




cos 9 Sin 9 
— sin 9 cos 9 




(4) 



where 9 denotes the Z — Z' mixing angle. In addition, the physical states of quarks could be 
related to the weak eigenstates by = Vi[l(^r)U'^ and -D^^^-j = Vdl{r)D^ in which U'^ = 

{u, c, if:), = {d,s,b), and the Vu^fi) and VDL(i?,) are the unitary matrices for diagonalizing 
weak states to physical eigenstates. The CKM matrix is defined by Vckm = VulVdl- -^^ ^ 
result, in terms of physical states the effective interactions for b ^ s decays could be written 
as 

Hz = ^sin^cos^ J2 E [BTefh,Px.sqYPx^ci]+h.c. , 

nz' = ^ cos^ 9 ' E [BTB-fb^,Px.s -qYPxA + h.c. (5) 

where Tiz and Tiz' express the effects oi Z — Z' mixing and Z' , respectively, the q could be 
M, (i, s and c quark, Xi=i,2 = L and i?, and 

= Tl - Qg sin^ 9w, ef = -Q, sin^ 9w, 

dX _ yD-XyD^ nx _yU~Xym (a) 
^DD — ^x D^x ' UU ~ ^x ^U^x ■ V / 

Here, the capital UU and DD in the subscript of parameter could be the flavors (u, c) 
and (d, s, b), respectively, and the 9^ is the Weinberg's angle. By current experimental 
data, it is known that the mixing angle 9 is limited to be less than O{10~^) ll(J|. If the 
mass of Z' is in the range of a few hundred GeV to 1 TeV, the dominant effects only come 
from the Z' exchange. Therefore, under this assumption, we will neglecte the contributions 
of Z — Z' mixing. 

According to the interactions of Eq. (0), the new effective Hamiltonian for b — sqq 
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decays could be written as 



H. 



z' 



cS 



Gp 
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VtlVts J2 { (i-'^v-A (aCs + AC. 



{qq)v-, 



+ ( AC5 + ACr^eg 



+ 



AC^ 



AC: 



,3 



{(iq)v+A 
{q<i)v-a 



+ {bs)v+A 



V+A 



(7) 



with 



AC 



3[5] 



AC3[5] 



AC 



9[7] 



ACg^y] 



2 f g^Mz 
'3 \giMz> 

2 f g^Mz 
'3 \giMz> 

4 f g^Mz 
"3 Ui^^z' 

4 /^2M£ 
"3 \giMz' 



V*Vt 



B^: [ B^jf + 2B 



tb^ts 



MR] 
'DD 



-B 



R* 
sb 



Buu +'^^DZ) ) 



1 



B^: ( 5 



B 



RIL]\ 
DD J 



(8) 



The expressions have been written as the four-fermion operators of the SM, shown in Eq. 

The operators associated with the new unprimed WCs ACs^sj^g are the same as SM. 
However, the operators associated with primed coefficients AC3 579 have different chirahty 
from those in the SM for 6 — s couphngs. That is, the flavor-changing (FC) Z' model provides 
different chiral flavor structures for FCNC processes. It has been found that the new effective 
WCs of Eq. (jSI) could be simplified if one assumes B^jj ~ —25^^ Although in general 
the assumption is unnecessary, for simplicity, we still impose the condition in our case. 
Hence, we get Ad'l ^ and 



AC 



9[7] 



ACgJy] 



4 



f 92Mz 
\9iMz' 

( 92MZ 
\giMz' 



TjL* r)L[R] 
^sb ^DD 1 



(9) 



Although the hadronic matrix elements, describing the B decaying to two final mesons 
through the effective Hamiltonian, depend on the chiral and color structures of four-fermion 
operators, we find that the associated effective WCs could be classified and reexpressed to 
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be more useful form by 



c 

NP 



IqNP 3 , iqj^p 3 ACg iqi^p 3 , ,g^p 3 ACy . 

with Cg^f^ = Cq(j) + AC9(7). The superscript g of Eq. (fTIH) denotes the corresponding flavor 
and Cq is its charge. Since the considering Z' model has the flavor structures which are the 
same as SM, to be more clear to understand the influence of new physics, we rewrite Eq. 
(imil to be 

q qSM 3 . „ q gSM , 3 ACg 

% = al + -CgACg, al = al + 

q qSM , 3 . „ q qSM , 3 AC7 , . 

«5 = «I + 2e'?^'^7, ^6 = «6 + 2''^1^- 

III. DECAY AMPLITUDES FOR B AND B ■k{p)K^*^ DECAYS 

To describe the amplitudes for B decays, we have to know not only the relevant effective 
weak interactions but also all possible topologies for the specific process. We display the 
general involving flavor diagrams for h —>■ sqq in Fig. ^ where (a) and (b) denote the 
emission topologies while (c) is the annihilation topology. The flavor q in Fig^a) and (b) 
is produced by gauge bosons and could be m, or d or s quark if the final states are the 
light mesons; however, q' stands for the spectator quark and could only be u or d quark, 
depending on the B meson being charged or neutral one. However, the role of q and q' in 
Fig.^c) is reversed so that g = u, or d or s is the spectator quark and g' = u or d is dictated 
by gauge interactions. We note that the presented flavor diagrams are based on the penguin 
operators of the SM. Except the different type of interactions at vertices, the flavor diagrams 
induced by new physics should be similar to those generated by the SM. In addition, since 
the matrix elements obtained by the Fierz transformation of 03^4 are the same as those of 
Oi 2i we don't further consider the matrix elements of tree operators. Hence, in terms of the 
effective interactions of the SM and those shown in Eq. ((Tj), the expressions of Eqs. fllO|) 
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\r V 

b , \ / , q b , Y . 



q q q 

(a) (b) 



(c) 

FIG. 1: The flavor diagrams for b sqq: (a) and (b) stand for the emission topologies while (c) 
is annihilation topology. 

and (fTT|) . and the flavor diagrams of Fig. ^ we can investigate the decay amplitudes for 

B ^ Kcj), B ^TiK, B ^ K*(l) and B pK*. 



A. Bd ^ i^V 

Although there are charged and neutral modes in i? ^ K(f) decays, since the differences 
of flavor diagrams in charged and neutral modes are only the parts of small tree annihilation, 
we only concentrate on the decay Bd K^cj). At quark level, the process is controlled by 
the decay h — > sss, therefore, q = s and q' = dm Fig. ^a) and (b), but they are reversed in 
Fig. ^c). Hence, the decay amplitude for B^ K^cj) is written as 

-Mi; = V:,Vu [UilF^K + fB{i2F^K4> + ^3^2^^)] (12) 

and 

4l - a + a , 42[3] - «4[6] + '^4[6] ) l-'-'^J 

with a* = + a| + and a'^^^ = a'^'^^ + a'^^^ + a'^^^ . The means the decay constant 
of meson and is deflned by (O|s7^s|0) = f^m^e^, in which and express the mass 
and polarization vector of (f) meson. The hadronic matrix element is from the diagrams 
(a) and (b). However, the F^^^ and F2j^^ come from the annihilation topology diagram (c). 
The detailed expressions of the hadronic matrix elements are given in the Appendix I VI (^1 
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The decay rate for B PV is written as 

^^G|m|^l^|2 (14) 
167r 

where Pc = m\m\{r'^ — l)/m\ with r = Pi ■ P2/ {''^1^2) is the momentum of the outgoing 
vector mesons. By the decay width, we can also define the direct CPA to be 

AcP = ^ (15) 
where f is the decay rate of antiparticle. 



B. B ^ttK 



There are four specific modes in i? ^ -nK decays. Since all BRs and CPA of n^K^ 
are observed well in experiments, we have to analyze all of them in detail. We begin the 
analysis from the decay B^ —>■ n^K^. According to the flavor diagrams Fig. ^a) and (c), 
the emission and annihilation topologies of the decay are described hj q = d and q' = u. 
Hence, the decay amplitude for B^ n^K^ decay could be expressed by 

where fx is the decay constant of kaon and defined by (0|s7^75(i|-ft'°) = fKP^J., = m^/rriB 
with m^^ being associated with (0|s75(i|K°) = fx'm'j^, and 

Ct = al- af^, CI = at- a'e^^- (17) 

The hadronic matrix elements Ff^, F^^j^ and F2^j^ are similar to those for B^j K^cf) and 
the detailed expressions are given in Appendix lVIBI In addition, we have a new contribution 
F^-kK which arises from the emission topologies of 06,8- 

Similar to B^ n^K^, we can obtain the decay amplitude of Bd 7i^K~^ easily by 
using g(g') = u{d) instead of g(g') = d{u). Thus, the decay amplitude for B^, n^K^ 
decay is written as 

-V:,VusfKaiFl^, (18) 

where we have included the tree contributions. As mentioned before, except the CKM 
matrix elements and effective WCs, the hadronic effects of tree are the same as the penguin 
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operators 03,4. Hence, the hadronic effects encountered in Bd 7[~K~^ decay are the same 
as those in n^K^. 

Next, we analyze the situation of B^ n'^K^ decay, one can easily find that besides 
the involving flavor diagrams appeared in the decays B^ tt^K^ and B^i n^K^, the 
diagram Fig. E^b) corresponding to the electroweak penguin contributions of the SM should 
be also included. Taking proper flavors for q and q', the decay amplitude for Bd n^K^ 
decay is given by 

-V:,V^sf.a,F^K, (19) 

with ( = - ai + - - {af^^ - a!^^^ + af^^ - a'^^^). We note that the new 
term fnFij^, corresponding to Fig. ^b), has opposite in sign to other terms. The reason 
comes from the flavor wave function of 7r° being {uu — dd) / \p2. The Fig. ^b) picks both 
components while others only take dd component. Since the tree contributions are color 
suppressed, the corresponding WC is ai- 

After introducing the decay amplitudes of B^ —>■ tt^K'^, Bd —>■ ir^K^ and Bd n'^K^, 
the amplitude for B^ 7T^K~^ decay could be immediately obtained as 

V2M'^0K. = V*V,s [MC^F^^ + QtkF^^) + ACF^V + fsiQFU + QF^)] 

-V:,VUfKaiFl + Ua^Fl^). (20) 

Clearly, the amplitudes shown in the first three decay modes all appear in the decay B^ — >■ 
n^K^ . That is, once one determines the first three decays, the decay B^ -k^K^ is also 
fixed. 



Bd^K 



For the production of two vector mesons in B decays, since both vector mesons carry 
spin degrees of freedom, the decay amplitudes are related to not only the longitudinal parts 



but also transverse parts. In terms of the notation of Ref. 
be expressed by 



2^, the amplitude M^^^ could 



2p 



(21) 
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with the convention e°^^^ = 1, where the superscript h is the hehcity, the subscript L 
stands for h = component while N and T express another two h = ±1 components, and 



el(t) ■ 62 (t) = 1 with t = ±1. Hence, each hehcity amphtude could be written as 



22 



Ho = mlMi 



H± = m%MN =F mv^mv^\/r'^ - IMt ■ (22) 
In addition, we can also write the amplitudes in terms of polarizations as 

Al = Ho A||(^) = i=(/7_±/7^). (23) 
Accordingly, the PFs can be defined as 

= I. 12 , , (^ = A|U), (24) 



Consequently, the decay rate for B ^ V2V1 is given by 

where = \Piz\ = \P2z\ is the momentum of either of the outgoing vector mesons. 

From Fig. Q it is easy to see that the associated flavor diagrams for B —>■ K*(j) are the 
same as those for B Kef) decays. Furthermore, as the results for the neutral and charged 
modes are expected to be similar by neglecting the small annihilation contributions from 
tree operators 0^2 appearing in the charged mode, which will be discussed in Sec. IIV Al We 
will concentrate on the neutral B decay. In terms of the distribution amplitudes of vector 
mesons, defined in Appendix IVI Al the decay amplitudes with various helicities defined by 
Eq. (PT|) are given by 



(26) 



where H = L, N, T and 



— U|| — a — a , — a + a 

42[3]L - 42[3]|| - a4[g] - a4jg] , 42[3]± - a^Q] + 04(6] • l^' j 

The definitions of a** and a'^^^ are the same as those for Bd — > K'^ip decay. The explicit 
expressions for could be referred to Appendix IVI Dl 
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D. B ^ pK* 

Since the quark compositions of p and K* mesons are the same as those of vr and K, 
respectively, the flavor diagrams for B —>■ pK* and B nK decays should be the same. 
However, due to {V{p)\qq'\0) oc ey(p) ■ p = in which the scalar vertex is arisen from the 
Fierz transformation of {V — A) ^ {V + A), the emitted factorizable contributions of four- 
fermion operators 06,8 ^-re vanished, i.e. have no contributions. Consequently, it could 
be expected that BRs of 5 — > pK* are smaller than those of B ^ ttK in the SM. 

Although there are four possible modes in the B —>■ pK* decays, we only concentrate 
on the decays B^ — > p'^K*^ and B^ p^K*^ that they have larger BRs. Following the 
definition of Eq. (j^H), the various helicity amplitudes for B^ p^K*'^ decay would be 
written as 

■^f+K*OH = Vtb^ts [fK*CtHFpH + IbCihPipK'H + fBQHF2pK*H\ ■ (28) 

The associated effective WCs are given by 

—1^1 — f'^ — /-I 

^1H=L,\\ ^2H=L,\\ ''1' ^3H=L,\\ 

/■g _ A? _ „9 I 'gNP _ q , 'qNP /qq\ 

'^1± — '^2± - O'A + O'i 5 S3_L - '^6 + ^6 • l^^J 

Similarly, the decay amplitude for B^ —>■ p~K*^ decay is written as 

■^f-K'+H = VtlVts [fK*ClH^pH + fBC2HPlpK*H + fBCtHP2pK*H\ 

-V:,VusajK^F;^. (30) 
The definitions of Cf(2) same as those for B — > ttK, expressed by Eq. (fTTj) . 

IV. NUMERICAL ANALYSIS 
A. Theoretical inputs 

To obtain numerical estimations, the values of theoretical parameters in the SM related 
to the weak interactions are taken as follows: Gp = 1.166 x 10^^ GeV^^, Vus = 0.224, 
Vts = 0.041, Vub = 3.5 X 10~^e~*'^^ with 03 = 72^. The decay constants of mesons are set 
to be U = 130, fx = 160, fn = 190, fP = 237(170), and = fP = 200(160) MeVs. 
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The lifetimes of charged and neutral B mesons are chosen as tb+ = 1-67 x 10~^^ s and 
tb^ = 1.56 X 10^^^ s, respectively. Since we use PQCD approach to calculate the hadronic 
matrix elements, we set the scale of weak WCs at /i = \/AmB ~ 1.5 GeV, therefore, the 
values of SM shown in Eq. (finjl are estimated to be 

ai = 1.07, aa = -0.028, 
^uSM ^ _o.o02, af *^ = -0.036, af^ = -0.008, af^^ = -0.055, 
af = 0.012, af = -0.036, af = -0.008, af = -0.056, (31) 

and In addition, in Table |l] we present the values of the hadronic effects, which 

are displayed in Sees. IIII AI - ITlI Dl and calculated by PQCD approach. From the table, we 



TABLE I: The values of factorizable amplitudes. 





771a 


rpa 

^2K<f> 




^2tv 


0.37 


(-9.88 + z7.54)10-^ 


-0.047 + i0.14 


0.24 


0.50 




rpa 


FhL 




Ffc*^ 


(0.39 + ^8.16)10-'' 


(1.99 -i3.36)10-2 


0.36 


0.06 


0.11 


rpa 


rpa 


rpa 


rpa 

^2K*4>L 


77a 

^2K''4,\\ 


(-1.4-^1.0)10-3 


(6.6 + ^6.5)10-4 


(-1.2-^6.4)10-3 


-0.03 + ^0.14 


0.03 -i0.02 










rpa 

^lpK*L 


0.06-fO.ll 


0.31 


0.04 


0.08 


(-2.3 - i5.4)10-3 


rpa 

^lpK''\\ 


rpa 

^lpK*± 


rpa 

^2pK*L 


rpa 

^2pK*\\ 


77a 

^2pK*l. 


(3.1 + i0.9)10-^ 


(-1.9 + ^2.9)10-3 


0.03 + i0.16 


(0.65 -i8.3)10-2 


0.01 - i0.17 



clearly see that the annihilation contributions from {V — A) ® {V — A) operators which 
correspond to F^pp and F^yy are negligible. To be more clear understanding the results, 
we use B PP decays to illustrate the property. For B PP decays, the factorized 
amplitude associated with the {V — A) ^ (V — A) interaction for annihilated topology can 
be expressed as 

{PiP2\qiYil - 75)^2 537^1 - l6)b\B), = -tfBiml - ml)Fl^^''\ml) (32) 
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where mi(2) are the masses of outgoing particles and -FQ^^^^(m^) corresponds to the time-hke 
form factor, defined by 



{0\qYl,b\BipB)) = ifBP% 
{Pi{pi)P2{p2)\qilM^) -- 



respectively, with q = pi — P2 and Q = pi + P2- From Eq. (jHSl), it is clear that if mi = m2, 
the factorized effects of annihilation topology vanish. However, the cancelation factor will 
be removed when the interactions correspond to {V — A) ® {V + A) operators 



B. Experimental inputs and predictions of the SM 

As mentioned before, the accuracies of some experimental data on BRs and CPAs are 
quite well, thus we could utilize these observed values to constrain the new parameters of 
the Z' model. To be more clear to know what the experimental inputs and the predictions 
are, in the following we definitely display the ranges of current experimental data for the 



inputs. Hence, taking the world averages with 2cr errors presented in Ref. j24|, the inputs 
of BRs are -B^ K^(t>i B —>■ K*^(j) and all B ttK decays, and their limits are taken to be 

7.3 < BR{Bd is:0)lO^ < 9.5, 8.6 < BR{Bd K*(p)10^ < 10.4, 

21.5 < BR{B^ 7r^K)m'^ < 26.7, 16.6 < BR{Bd tt^K^)10^ < 19.8 

9.5 < BR{Bd^ 7r°K)10^ < 13.5, W.5 < BR{B^ ^ 7r°K^)10^ < 13.7. (34) 

Moreover, we also take into account the ratios of BRs, defined by 

BRjBd ^ n^K^) ^ 2BR{B^ n^K^) ^ BR{Ba ^ vr^ir^) 

' TB, BR{B^ -> 7r±K) ' ' BR{B^ ^ n^K) ' " 2BR{Bd ir^K) 



as 0.76 < Ri < 0.88, 0.91 < Rc < 1.09 and 0.74 < R„ < 0.88 [11<]. Since there are no 
measurements on the CPAs oi B ^ K^*'^(f), we artificially set the limits as < \Acp{B — >■ 
K^*^(l))\ < 0.05 in which the CPAs vanish in the SM. Other limits from data are taken as 
< \Acp{B+ 7i+K^)\ < 5%, 7.1% < \Acp{Bd 7t-K+)\ < 14.7%, and < \Acp{B+ 
TT^K^)\ < 10%. Because there is no any significant information on the CPA of B^ tt^K^, 
we leave the value as our prediction. In addition, we also take the longitudinal polarization 
Rl of Bd K*°(t) as the input and the limit is chosen to be 44% < RiiB^ K*°(t)) < 57%. 
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Before going on discussing the contributions of the Z' model, it is worth knowing the SM 

results which are based on the taken values in Sec. IIV Al Hence, the SM predictions on BRs 
are 

BR(Bd. K^) = 8.98 x 10"^ BR{Bd K*(j)) = 12.9 x 10~^ 

BR{B^ p^K*) = 15.3 X 10"^ BR{Bd p^K*^) = 13.4 x 10-^ 

BR{B^ n^K) = 22.2 x 10"^ BR{Bd n^K^) = 19.0 x 10"^ 

BR{Bd n^K) = 7.87 x 10"^ BR{B^ tt^K^) = 12.5 x 10"^; (36) 

the ratios of BRs are estimated by Rf^' = 0.92, R^^^ = 1.18 and R^^ = 1.22; and the 
predictions on CPAs are 

AcpiBd ^ p^/r±) = 22.0%, AcpiBd ^ n^K^) = -12.1%, 

AcpiBd ^ n'^K) = -1.35%, Acp(5± ^ n'^K^) = -8.3%, (37) 

where the vanished CPAs are not shown. Moreover, the estimations of the various PFs for 
VV modes are also given to be 

RL{Bd K*^<P) = 0.71, R^{Bd ^ = 0.15, 

Rl{B+ ^ p+ir*°) = 0.72, R^{B+ p+K*^) = 0.13, 

RiiBd ^ p-K*+) = 0.52, R^{Bd ^ p'K*+) = 0.22. (38) 

According to our estimations, we see that compared to the data, B^ K*(f){Bd -n^K) 
has larger (smaller) BR, the ratios Ri^c,n don't fit the data well, and Rl of Bd K*^(f) is 
much larger than observations. We also find RlIB^ p'K*^) could be around 50%. For 
displaying the influence of different scales, in Fig. |21 we present the correlations between BRs 
in K(l) and K*(j) modes, Rl,i. and BR{Bd K*(f)), and Rl and R±, where the circle, square, 
diamond and triangle-up symbols stand for the results of /i = 1.3, 1.5, 2.0 and 4.0 GeVs, 
respectively. The error bars presented in the figures are the world averages with 2a errors. 
Similarly, we also show the SM predictions on the CPAs of -B — > ttK and the corresponding 
BRs in Fig. El 



C. Results of the Z' model on B ^ (t'K* , B ttK and B pK* decays 



Before performing the numerical calculations, we first discuss the allowed regions of new 
effects which are from ACgj-,j . According to the results of Refs. 3, ll^ , it is known that 
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1.5, 2.0 and 4.0 GeVs, respectively. The error bars are the world averages with 2a errors. 
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the unknown parameters, defined by = {g2mz/ gimz'YB^^Bf^/V*fjVts with X = L, R, 
have been hmited to be |^^"^| < 0.02 in which m'^ is at TeV-scale. That is, if we assume 
that B^fj ~ B^f^ and -B^^' ~ -^dd'' consequently we obtain |ACgj^j| < 0.08. In the following 
analysis, we will take this value as the upper bound of new effects. Since ACgj^^j in general 
are complex, we totally have eight parameters for each b — * sqq with q = {u, d) and 
h — * sss decays. In principle, the eight free parameters for h sqq could be fixed by the 
eight chosen measurements such as four BRs and four CPAs in i? ^ tiK decays. Then, 
using the constrained parameters we can make predictions on B ^ pK*. Although there 
are no eight measurements related to 6 — > sss directly, however, due to the new effects in 
Eq. (fT^ for B Kcj) being different from that in Eq. (j27j) for ^i(2)l of i? ^ K*(j), we find 
that when the data of Kcj) and K*(j) are considered simultaneously, the unknowns has been 
strictly constrained. For convenience, we parameterize the unknowns to be ACg = r^LLC*"^^^, 
AC7 = riLRc'^^^, AC^ = riRLc'^^^, AC'^ = tirrc^'^R^ so that |r/xy| < 0.08 and < 0xy < 27r 
with X and Y each being L or R. 

Now, we could investigate the contributions of the Z' model to the considering processes. 
At first, we study the decays governed by 6 — sss. It has been known that in terms of fiavor 
diagrams of Fig. ^ all effects contributing to B^ —>■ Kcj) will also influence on B^ K*(j>. 
It could be expected that in SM-like models, to reduce the longitudinal polarization Rl 
of Bii K*(j) will also lower the BR of i? ^ Kcj). We flnd that based on the hadronic 
values of Table HI if we tune ACgj^j = 0, there are no solutions for the ACgir] to satisfy 
the data of BR{Bd K^*^cj)) and RL{Bd K*cj)) at the same time. And also, if we set 
AC9[7] = ACgj7j or AC9[7] = AC^^gj etc, no possible solutions are found. That is, in order 
to flt the current experimental data, f]xY{4>XY) cannot have simple relationship for different 
X and Y. Hence, by taking each 7]xy < 0.08 and each cpxv = [— tt, tt] and including the 
limits of Eq. (jH^ and 44% < Rl^B^ K*^cj)) < 57%, we present the possible solutions in 
Fig. El By Fig. El^a), we could see the correlation of BR between Kcj) and K*cj). From the 
Fig. m^b) and (c), we see clearly how the changes of Rl{i) are associated with the BR of 
K*cj). We also present the correlation of Rl and -R_l in Fig. |31^d). According to these results, 
it could be concluded that Z' model which provides the left- and right-handed couplings 
could solve the anomalies of small Rl{B — > K*cj>). In addition, the Z' model also provides 
the room for large R_l{B K*cj>), say above 25%, in which R± of the SM is around 16%. 
As comparisons, we also show the results of fi = 1.3 GeV in Fig. El Clearly, more solutions 
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are allowed. We note that no solution can be found when /x > 1.5 GeV. 
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Concerning the processes dictated by the decays b — ^ sqq, we also find that if we tune 

ACyjg] etc, no possible solutions for simultane- 



AC^j7] = 0, or AC9[7] 



AC^j7], or AC9[7] 
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ously matching the data are found. Hence, all unknowns should be regarded as independent 
of parameters. Following the formulas introduced in Sec. IIIIBl the constraints of Eqs. (jM)) 
and (j35|) as well as the bounds of CPA, we display the results in Fig. [Ul From the figure, 
we could see that the CPA of tt^K could be as large as 4% while it vanishes in the 

SM. To satisfy all current experimental data, the CPA of tt^K should be smaller than 

— 10%. Since the CPA of tt^K^ has a good accurate measurement, if one could further 



confirm that the magnitude of CPA for B^ tt^K is small, say less than 



we could 



conclude that the large CPA of B^ — > n^K could be a very good evidence to display the 
existence of new physics, where the SM prediction is only around —3%. Furthermore, by 
Fig. Efd), we also see that the CPA of B"^ — > tt^K^ could be much smaller than that of 
B^ TT^K"^, in which they should have similar values in the SM. And also, the results show 
that the CPAs of tt^K^ and tt^K^ favor to be opposite in sign but the SM predicts the same 
sign. As mentioned in the end of Sec. IIIIBl when the decay amplitudes for B^ n^K, 
Bd n^K^ and Bd n^K decays are determined, those for B^ n^K'^ decays are also 
fixed. Therefore, the sign difference could be also as the clear evidence that new physics 
exists. In Fig. [3 we also presented the results with = 1.3 GeV. Since the data of B ^ tiK 
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have better accuracies, we find that no possible solution appears when the scale is smaller 
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Finally, we discuss the contributions of the Z' model to the decays B pK*. By 
the analysis in Sec. IIIIDl we know that except the transverse parts, the weak WCs for 
longitudinal polarization Rl of B —>■ pK* should be the same for the decays B ttK, 
i.e. they have the same weak effective WCs Ci,2, as shown in Eqs. (fTTj) and (jSH))- This 
is because the final states in both processes have the same parity properties. However, 
the case encountered in the decays B K^*'>(j) is different because the parity properties 
of final state Kcj) are different. Hence, the values constrained by i? ^ -kK could directly 
make predictions on B ^ pK*. We present the results of B^ p^K* and B^ p'^K*^ 
in Figs. IHl and 01 respectively. Since the observed BR of B^ tx^K has reached a good 
accuracy, in FiglHfa) and FigEfa) we show how the BRs are associated with BR(i?^ — >■ 
tt'^K). Moreover, we display the CPA, Rl and R± versus the corresponding BR in (b), 
(c) and (d) diagrams of both figures, respectively. We note that the observed BR(_Ri) of 
5^ p^K* by BABAR and BELLE are not consistent each other. The former observes 
.7.0 ± 2.9 ± 2.0(0.79 ± 0.08 ± 0.04) while the latter is 8.9 ± 1.7 ± L0(0.43 ± 0.11^°:°^) 
2fil |. By the Fig. |Hl we could see clearly that (1) 5^ p^K* can have sizable CPA in which 
it vanishes in the SM; (2) Rl could be less than 0.60 while the corresponding BR is above 
15 X 10^^; (3) the solutions of small R±_ exist, i.e. R\\ » R± where the prediction of SM 
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is R\\ ~ R±. As for the results of Bd p^K*^ shown in Fig. 13 due to just hke the case of 
Bd Ti^K^ which the results with new effects are similar to the SM, we expect that the 
derivations from the SM are not too much. Hence, we could summarize the favorable ranges 
of BRs, CPAs, Rl, andi?x for (5± p^K*, B^ p^K*^) are (17.1±3.9, 10.0±2.0) x 10"^ 
(3 ± 5, 21 ± 7)%, (0.66 ± 0.10, 0.44 ± 0.08) and (0.14 ± 0.10, 0.25 ± 0.09), respectively. 



24 

I 

o 

23 

T 

+1 

CQ 

^ 22 
CQ 



12 



1 1 


1 1 1 1 1 1 


1 1 1 


1 1 1 


















1 1 


1 1 1 1 1 1 


, 1 , 


, 1 ,- 



15 



18 



21 



BR(B"^ p"K )10 



-6 




BR(B-^p-K )10 



15 [ 

10 ; 



5 \ 

41 

« 

u 

< -5 



12 



0.4 r 

"Q 0.3 ; 

+1 

Q. 

T 0.2: 

41 



0.1 



12 



itOliiliSi:;;; 



15 



18 



21 



BR(B"^p"K )10 



(d): 



15 



18 



21 



BR(B"^p"K )10 



FIG. 8: (a) correlation of BR (in units of 10 ^) between p^K* and n^K; (b), (c) and (d) denote 
the correlations between (CPA, i?_|_) and the BR, respectively. 



V. SUMMARY 

We have studied the effects of nonuniversal Z' model on the processes dictated by the 
h sqq decays with q = u, d, and s. By using the PQCD approach, we calculate the needed 
hadronic matrix elements. For B K^*^(f) decays, we find that their BRs and the Rl of 
Bd K*(f) have provided strict constraints on the new parameters. After marching the 
currents data, we find the -R_l of Bd — > K*^(f) favors to be larger than 25%. For B ttK 
decays, by requiring that the magnitude of Acp{B^ tt^K) is less than 5% and all BRs 
satisfy the current observations, we find that the magnitude of CPA of Bd tt^K should 
be larger than 10% but sign is the same as SM. Meanwhile, the CPA of B^ tt^K^ could 
be as low as few percent which is indicated by the current experiments. Moreover, we also 
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FIG. 9: The legend is the same as Fig. IHlbut for Bd p^K*^. 

obtain that the CPA of —>■ -k^K^ is opposite in sign to the SM. 

In sum, to satisfy current data, the new left- and right-handed couplings have to be 
included simultaneously. It is clear that the FC Z' model provides the needed couplings 
naturally. With more physical observations and accurate data by B factories, we could 
further examine the effects of nonuniversal Z' model. 
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VI. APPENDIX: DISTRIBUTION AMPLITUDES AND DECAY AMPLITUDES 



A. Distribution amplitudes 

We describe the spin structures of meson to be 



{P{p)\q2p{z)qUm) = -^== / rfa;e"^-^{75i'$p(a:)+75<$?,(x) 

V ^^Vc Jo 

+m°,75(/^+ ^--IWp]^p (39) 
for pseudoscalar, where p = (p"*", 0, 0_l), = (1, 0, Oj_), and n_ = (0, 1, 0_l); and 

1 

\'ZI\c Jo 

+ ^^zeT^,,.757^eV^"*y (a:) | (40) 

(T) 

for vector meson. The notations $ p and denote the twist-2 wave functions while 
^p*^ and (l)^'^'"'" stand for the twist-3 wave functions of pseudoscalar and vector meson, 
respectively. Their explicit expressions could be found in Refs. 



B. Hard functions for B ^ P2P1 decays 

In terms of the spin structures of mesons defined by Appendix IVI Al we write the 
factorizable amplitudes ioi B ^ P transition form factors and B PP annihilations as 

F^p = S-n-CpMp / dxidxs / bidbib-sdb-s^B {xi,bi) 
Jo Jo 

X { [(1 + X3) <^P (X3) + rp(l - 2x3) (0:3) + n (X3))] (t«) 

X he{xi,X3,bi,b3) +2rp<l>''p{xs)Ee{t';^^)K{xs,Xi,b3,bi)} , (41) 

F2P = IQnCpM'^rp / dxidx^ / bidbib^dbs^^p ixi,bi) 
Jo Jo 

X { [<l>p (X3) + rp ((2 + xsWp (X3) - 0:2$^ (X3))] (t«) 

X /le(Xi,a;3,6i,63) +2rp<l>^(x3)Ee (4^)) /le(x3,Xi,63,6i)} , (42) 
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n 1 noo 

= -^T^CpMl / dx2dx^ / h2dh2hdh 
Jo Jo 

X { [s3$P,(s2)$P,(l - X3) + 2rp,rp,<l'^^(x2) ((1 + a;3)$?,^(l - X3) 

+ (1 - X3)$^,(l - X3))] (4^0 /ia(x2, X3, &2, ^s) 

- [X2$P,(X2)$P2(1 -X3) + 2rp,rp,<l>^^(l -X3) ((1 + X2)$?,^(x2) 

-(1 -X2)$?,,(X2))] ^„ (ti')) /la(x3,X2,63,&2)} , (43) 



/•l POO 

^2P2Pi = 167rC^M| / dx2dx3 / b2db2hdb3 
Jo Jo 

X { [rp,a;3$P,(a;2) - 0:3) + $^,(1 - 0:3)) + 2rp,^^p^{x2)^P,{l - X3) 

xEa (ti^^) ha{x2,X3,b2,b3) + [x2rp, ($?,^(x2) - $p,(a;2)) $P2(1 - X3) 
+2rp,$p,(x2)$^,(l -X3)] Ea (ti')) K{x3,X2,b3,b2)} (44) 

with iriPi^^/mB-i where the hard functions he{a) are given by 

K{xi,X3,biM) = Ko (y/x^niBbi) Stix3) 

X [9{bi - 63)^0 i^/x^mBbl) Jq (Va^^is^s) 

+61(63 - ^1)^0 (v^^^B^s) ^ (y/x^ruBbi)] , (45) 

ha{x2,X3,b2,b3) = |^y^ ^^(S^^ {y/ X2X3mBb2) 5't(x3) 

X 6'(62 - 63)-f^o^'' {Vx^rnBb2) Jo (v^^^b^s) 

+^(63 - 62)i^i'^v^mB63) -^0 (v^mp62)] • (46) 

The evolution factor -E'e(a) are defined as 

(t) = as (t) Sb (t) Sp {t) , 
E,{t) = as{t)Sp,{t)Sp,{t) (47) 

where SmU) denote the Sudakov factor of M-meson, the exphcit expressions could be found 
in Ref. l22l and the references therein. 
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C. Hard functions for B PV decays 



Similarly, the factorizable amplitudes for B — > PV modes are given to be 

/» 1 /*oo 

F^py = StiCfMI / dx2dx3 / b2db2hdh 
Jo Jo 

X {[X3^v(x2)^p(l - X3) + 2rprv^'y{x2) ((1 + X3)^>(1 - X3) 

+ (1 - X3Wp{l - X3))] Ea{t^^^)ha{x2, Xg, 62, ^s) " [X2^v{x2)^ p{l 
+2TPTV ((1 + X2WV{X2) - (1 - X2)^'y{,X2)) - X3)] 

xK(ii^^)/ia(a;3,a:2,&3,M| > 
io io 

X {[rpX3$y(a;2) - X3) + $^(1 - X3)) + 2ry$t'(^2)*p(l - 

X£;„(i(l))/i„(x2, X3, b2, 63) + [2rp$y(x2)$?,(l - X3) 

+X2rv {^v{^2) - ^vi^2)) $P(1 - X3)] Eaitf^)haix3,X2,b3,b2)} 

with rv — mv/mB- 

D. Hard functions for B V2V1 decays 

The needed factorizable amplitudes for VV modes are given by 

/■I poo 

F^^ = SnCpMl / dxidx3 / bidbib3db3^B{xi,bi) 
Jo Jo 

X { [(1 + 0:3)$^, (0:3) + rv,{l - 2X3){^U{X3) + ^'v,i^3))] 
XE%t^P)K{x^,X3MM) 

+2rv,^\r^{x3)E\tf)K{x3,x^,b3M)] , 

pi poo 

Fy^N = SnCpMl / dxidxs / bidbib3db3^B{xi,bi) 
Jo Jo 

xrv, {[^iMi) + 2rv,<^l^{x3)+rvM'^ki^-i) ^U^s))] 

XE\t^^^)he{xi,X3MM) 

+rvAn,{^^) + n,{^^W{tf'^)he{x3,x,MM)] , 
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/•l POO 

Fy^j, = WttCfMI / dxidx^ / bidbib3db3^B{xi,bi) 
Jo Jo 

x£;^(4i))/ie(xi,X3,6i,63) 

+rvA^U^3) + nAx3Wit^e^)he{x3,Xl,b3,bi)} , (52) 
p1 poo 

= StcCfMI / dx2dx3 / b2db2b3db3 
Jo Jo 

X {[-X3$yi(x2)$y2(l - X3) + 2rv,rv2^v,ix2){{l - X3)^U{1 - X3) 

+ (1 + X3)^'y,(l - X3))] E-(tl^^)ha(x2, X3, 62, ^s) 

+ [x2$yi(a;2)$Mi(l - ^3) + 2rv,rv2<l>v^{l - X3){{1 - 2:2)$^^ (2:2) 

-{l+X2)^'y^{x2))]E^{tf^)ha{x3,X2,b3,b2)} , (53) 
pi poo 

^iVaViiv = -SnCpMl / dx2dx3 / 62^6263^63 

Jo Jo 

XTv.Tv, { [(1 + X3){n,{x2)n,{^ - X3) + $^,(X2)$^,(1 - X3)) 

+(1 -a;3)($^y,(a:2)$^,(l - X3) + $^^(a;2)$^,(l - 0:3))] i?'^(tii))/i.(a;2, 3:3, 62, 63) 
- [(1 + X2){nS^2)n,iy - ^3) + n,{^2)n,{^ - X3)) 

-(1 - X2Wv,{x2Wv,{^ - X3) + n,{^2Wv,{^ - X3))\ 

XE%t!^^)ha{x3,X2,b3,b2)} , (54) 

/•I /-oo 
^iVaViT = -IQttCfMI / dx2dX3 / 62^6263(^63 

Jo Jo 

xry.rv., { [(1 - X3){^l,{x2Wv,{l - X3) + (^2)^^,(1 - Xg)) 

+ (1 + X3)m,ix2)<^U^ - X3) + $^^(X2)$^,(1 - X3))] E%t'^^^)ha{x2, X3, 62, 63) 

+ [(1 - X2){<^"y^ix2)<^U^ - X3) + X3)) 
-(1 + X2){n,{x2)<^V,{l - X3) + ^t^M2WvAl - X3))\ 

XE-{t^^^)ha(x3,X2,b3,b2)} , (55) 

pi poo 

F2v,v,L = IGnCpMl / dx2dx3 / b2db2b3db3 

Jo Jo 

X { [rv2X3$yi(x2)($y,(l - X3) + $f.^(l - X3)) - 2ry,$f.^(x2)$y2(l " 2^3)] 

xE'^(tW)/i„(x2,a;3,62,63) 

+ [ry,X2($yj(x2) - $^^(x2))$y2(l - a;3) + 2ry,$yi(x2)$y,(l - 0:3)] 

x£;"(i(2))/i„(x3, 0:2,63,62)} , (56) 
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p1 poo 

^mViN = 167rC^M| / dx2dx3 / h2dh2hdh-i 

Jo Jo 



{ry,($^y^(x2) + $^,(X2))$^,(1 - X3)£;"(ti^))/la(x2,X3, 62,63) 
+ ry,$^^(x2)($i:.,(l-X3)-$^,(l-X3))i?^(t«)/la(x3,X2,63,62)} , (57) 



p i poo 

Fw.v.T = 327rCpM| / ^0:2^X3 / 62^6263^63 

^0 Jo 

X {ry,($;:.^(x2)+$^,(x2))$^,(l-X3)^''(ti^))/la(x2,X3,62,63) 

+ ry.$y,(x2)($y,(l-X3)-$^,(l-X3))i?''(ti'))/la(x3,X2,63,62)} . (58) 

We define ry. = mvJmB- 
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